Coherent quantum control of A-atoms through the stochastic hmit 
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We investigate, using the stochastic limit method, the coherent quantum control of a 3-level atom 
in A-configuration interacting with two laser fields. We prove that, in the generic situation, this 
interaction entangles the two lower energy levels of the atom into a single qubit, i.e. it drives at 
an exponentially fast rate the atom to a stationary state which is a coherent superposition of the 
two lower levels. By applying to the atom two laser fields with appropriately chosen intensities, 
one can create, in principle, any superposition of the two levels. Thus relaxation is not necessarily 
synonymous of decoherence. 

PACS numbers: 32.80.Qk, 42.50.Lc 



I. INTRODUCTION 

Preparation of atoms and molecules in a predefined 
state plays an important role in modern atomic and 
molecular physics, in particular in atom optics and quan- 
tum information. The difficulty of the problem consists 
in the fact that, in order to control a system, one has 
to interact with it, but interaction introduces dissipative 
effects and hence, at least generically, decoherence. In 
this note we prove, in specific but important and phys- 
ically realizable example, that in some cases dissipation 
can generate coherence. 

One of the ways to drive the system, atom or molecule, 
to the desired state is to exploit its interaction with laser 
pulses, i.e., to use the coherent laser control or laser- 
induced population transfer P, In this approach 
monochromatized and near-resonant with the atomic 
Bohr frequencies radiation fields are used to force the 
system to the final state. The laser coherent control tech- 
niques can be used in applications to laser cooling based 
on coherent population trapping j^j , quantum computing 
with trapped ions 0,0, etc. 

The stochastic limit method ^ was applied in ^ to 
study the phenomenon of coherent population trapping 
in an atom in lambda-configuration with a doubly degen- 
erate ground state. It was found that the action of the 
field drives the atom to a 1-parameter family of station- 
ary states so that the choice of one state in this family is 
uniquely determined by the initial state of the atom and 
by the initial state of the field. 

In T| the initial state of the field was chosen to be 
an arbitrary mean zero, gauge invariant Gaussian state, 
e.g. the Fock vacuum or an equilibrium state at any tem- 
perature. On the other hand, in the usual experiments 
on coherent population trapping, the atom is driven by 
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two laser beams, i.e. coherent states, resonant with the 
atomic frequencies and, up to now, there is no evidence 
of this 1-parameter family of invariant states 

In the present paper we prove that, for a three- level 
non-degenerate atom, the field drives the atom to a pure 
state, which is a superposition of the two lowest energy 
levels and that, applying to the atom two laser fields 
with appropriately chosen intensities one obtains a single 
superposition. In the case of a degenerate atom we show 
that, if one starts from a coherent state, then the family 
of invariant states is destroyed and the field drives the 
atom to the dark state. 

The results of the present paper, combined with those 
of also suggest that the emergence of the dark state 
could be experimentally realized not only by tuning two 
lasers, but also in a generic equilibrium state, by prepar- 
ing the state of the atom in the domain of attraction of 
the dark state. 

The present result is obtained as a consequence of a 
general property of the stochastic limit, which we find for 
the first time in the paper, according to which the effect 
of replacing a mean zero, gauge invariant state of the field 
by a coherent one, amounts to the addition of a hamil- 
tonian term to the master equation while the dissipative 
part of the generator remains the same [cf. H12l) - (|14|l ]. 

This general property is proved in sections 2-4 of the 
present paper for an arbitrary atom. In sections 5 and 
6 we specialize to the case of 2- and 3-level atoms, find 
for these cases the explicit form of the stationary atomic 
states [ (|17|) for a two-level, (|21ll and H22() for a three-level 
atom], and prove the results mentioned in the present 
introduction. The master equation we find in the 3-level 
case coincides with the optical Bloch equation consid- 
ered by Arimondo in Q in the absence of scattering 
[Eq. (2.7)] with the only difference that we obtain the 
equation and the explicit formula for the coefficients Qj 
[cf. equations © and (|20|) below] without assuming that 
the field is classical. Instead, we start from the micro- 
scopic quantum dynamics of the total atom+radiation 
system. 

The approach of the present paper to study the coher- 
ent quantum control or laser induced population trans- 
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fer is based on the derivation, using the stochastic hmit 
method, of a quantum white noise equation approxi- 
mating the total dynamics in the weak coupHng regime 
(Eq. ©, for a rigorous treatment see [^). The pro- 
gramme to exploit exponentially fast decoherence as a 
control tool to drive an atom into a pre-assigned state 
was formulated in 'lO', 'ill > where the control parameter 
is the interaction Hamiltonian. In the approach of the 
present paper the control parameter is the state of the 
radiation, which can be easily controlled in experiments 
and, instead of decoherence, it uses relaxation, due to dis- 
sipative dynamics, to create quantum coherence between 
the atomic states. 



II. AN ATOM IN A LASER FIELD 

The dynamics of an atom interacting with radiation 
is determined on the microscopic level by the interaction 
Hamiltonian and by the reference state of radiation. 

The Hamiltonian describing an atom interacting with 
radiation is given by the sum of the free and interaction 
terms 



int 



Ha ® 1 + 1 ® ffR + AiJi 



int ) 



where Hf^ and are free Hamiltonians of the atom 
and radiation, ffint interaction Hamiltonian, and A is the 
coupling constant. The free Hamiltonian of the atom has 
discrete spectrum 

Ha = ^£nPn, 
n 

where e„ is an eigenvalue and P„ is the corresponding 
projection. The free Hamiltonian of the radiation is 

H^^ j dkc^(k)a+(k)a(k), c^(k) = |k|, 

where the creation and annihilation operators a^(k), 
a(k) describe creation and annihilation of photons with 
momentum k. The interaction Hamiltonian has the form 

i?i„t = i{D (g> a+{g) - D+ (g) a{g)), 

where the operators D and describe the transitions 
between the atomic levels, a^{g) = / d'kg{\i)a^ (k) is the 
smeared creation operator, and the formfactor g(\i) de- 
scribes the coupling of the atom with the k mode of the 
field. 

The state of the radiation, which corresponds to a long 
time laser pulse of frequency uj and intensity / is a co- 
herent state, which is determined by the coherent vector 



where $o is the vacuum, W{-) the Weyl operator, which 
is defined as W{f) = exp[i(a(/) + a+(/))], St = e**"^') 



the one photon free evolution, and [S'/A^,T/A^] the time 
interval in which the laser pulse is active. This vector \ 
is an eigenvector of the annihilation operator: 



X- 



The eigenvalue 



ca = A 



J dt J dkg*(k)/(k)e 



S/A2 



is determined by the intensity of the laser / and the form 
factor g. The pulse starts at time 5/A^ and ends at 
time Tj}?. Since the coupling constant A is small the 
duration of the pulse, being of order A~^, is large. The 
function /(k) describes the amplitude of the input field 
with momentum k. 

The state, which corresponds to several laser pulses of 
frequencies wi, . . . , intensities /i, . . . , /„, and acting 
during time intervals [5*1, Ti], . . . , [S'„, T„], is determined 
by the coherent vector 



Ti/A" 



I A^ / Ste-'*'^^ fidt I $0. (1) 

1 = 1 c 



The initial state of the total system is supposed to be 
factorized 

UJ = uip ^ ujx = Tr (p •) (g) ujx, 

where Up is the state of the atom determined by a density 
matrix p and uJx is the state of the radiation determined 
by the coherent vector "ifx- ^\ = (^a, -^a}- 

The dynamics of the total system is determined by the 
evolution operator which in interaction picture is defined 
as 



Ux{t)^e 



where the index denotes the dependence on the coupling 
constant. The evolution operator determines the Heisen- 
berg evolution of any observable of the atom, say X, as 

Xxit) - U+{t)XUxit). 

The reduced dynamics of the observable is obtained by 
taking average over the state of the radiation 

{Xxit)) {^x,U+it)XUxmx) 

and, by duality, it can be equivalently described in terms 
of the reduced density matrix px (t) , defined by the fol- 
lowing equality: 



= Tr 



TvA{pxit)X) = TrA{p{Xxm 
rA {TVr [[/a (i) (/5 ® I^a) A l) C/a+ W] ^} , 
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where Tia denotes the trace over the atomic degrees of 
freedom and TrR the partial trace over the field degrees 
of freedom. 

The first step of the present paper is to give a non 
fenomenological derivation, from the first principles, of 
the quantum master equation for the reduced density 
matrix of the atom using the stochastic limit method. 
This is done in the following two sections. 



III. THE STOCHASTIC LIMIT 

It is impossible to find, for a general coupling D, the 
explicit form of the dynamics for the reduced density 
matrix for realistic models, while the weakness of the in- 
teraction suggests some approximations. For small time 
the dynamics can be effectively studied by using the per- 
turbation series. However, the n-th term of the series be- 
haves like (X^t)" and the approximation with the lowest 
order terms of the series becomes invalid for time which 
is large enough. 

Hence the study of the long time dynamics requires 
another approach. Such an approach to study the dy- 
namics on the long time scale (~ A~^) is the stochastic 
limit method 0|. In the stochastic limit one considers 
the long time dynamics (on the time scale of order A^^) 
for a system with weak interaction. Mathematically this 
means that one takes the limit as the coupling constant 
goes to zero, A — > 0, time goes to infinity, t — > cx), but 
in such a way that the quantity X^t remains fixed. In 
this limit the dynamics of the total system is given by 
the solution of a quantum white noise or stochastic dif- 
ferential equation (which is a unitary adapted process). 
Since the interaction of an atom with radiation is weak 
one can apply the stochastic limit procedure to study the 
long time dynamics of this system. 

On the long time scale the behavior of the exact re- 
duced density matrix is approximated by the limiting 
density matrix 

pW- limPA(t/A2) 

A — >u 

SO that p\{t) ~ p{X^t). The quantum master equation for 
the limiting density matrix in the case when the radiation 
is in a Gibbs state was derived and discussed in many pa- 
pers (see and references therein). The purpose of the 
present paper is to derive the quantum master equation 
for the limiting density matrix p{t) in the case the radi- 
ation is in a coherent state and to study this equation 
for particular, but important cases of two and three-level 
atoms. 

The evolution operator after the time rescaling satisfies 
the equation 



dUx{t/\^) 
dt 



where B is the set of all Bohr frequencies of the atom 
(spectrum of the free atomic Liouvillian i[HA, •], i.e. the 
set oi uj — Sn — £rrn whcrc Em Em are eigenvalues of the 
free atomic Hamiltonian), for each lu B 

ax At) = \j dk5*(k)e-'(-('^)~'^)*/^'a(k) (2) 

is the rescaled time dependent annihilation operator, and 
the operator 



E 



PmDPn 



(3) 



describes transitions between atom levels with energies 
En and Em with the energy difference En — £m = i^- 

In the stochastic limit one first proves that the time 
rescaled creation and annihilation operators converge, as 
A 0, to a quantum white noise (for details of the 
stochastic limit procedure see 6]): 

lini aA,c^(i) = b^{t), 

where the quantum white noise operators baj{t) are S- 
correlatcd in time and satisfy the commutation relations 



[bUt),b+{t')] - 6{t' - t)S^,^,2Rej^. (4) 
The complex number 

l5(k)P 



weB 



1^ = / — (5) 

J z((jj(k) — w — «0) 

is the generalized susceptibility and its real part 

Re7^ = 7r y"dk|g(k)|2,5(a;(k) -w) 

gives the decay rate of the w-transition (cf. 0, sect. 
4.20). The Kronecker 5-symbol in Q indicates the 
mutual independence of the white noise operators for dif- 
ferent Bohr frequencies. 

This allows us to derive the equation for the limiting 
evolution operator 

Ut := lim Uxit/X^). 
That is the white noise Schrodinger equation 



with the white noise Hamiltonian 

Vit) - i ^ (i?^ ® b+it) -D+® bUt) 
uieB 



Let X be any observable of the atom and Xt = 
U^XUt = limA^o ArA(t/A^) its limiting time evolution. 
Using the stochastic golden rule of the stochastic limit 
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(cf. 01, sect. 5.9), one gets the quantum Langevin equa- 
tion for Xt, that is a normally ordered white noise differ- 
ential equation: 



dXt 
dt 



U+Q{X)Ut 



+ 



h+{t)U+L+{X)Ut + U+L^{X)Uth^{t)) (7) 



(such an equation is equivalent to a certain quantum 
stochastic differential equation). Here 

Q{X) =J2{'^^^l^D+XD^-^^XD+D^-%D+D^X 

is the generator of a quantum Markov semigroup and, in 
the notation J^Jl, 



L+{X) = [X,D^ 



LUX)^[D+,X]. 



The normally ordered form of the quantum Langevin 
equation, when the annihilation white noise operators 
are put on the right of the evolution operator and cre- 
ation operators are put on the left, is very convenient for 
study the reduced dynamics in the case the radiation is 
in a coherent state. In order to get the reduced dynam- 
ics one just has to average equation (jT)) over the state of 
radiation. 

The coherent state of the radiation is approximated 
by the limiting coherent vector 

^ = lim vl/^. 

This vector is an eigenvector of the quantum white noise 
annihilation operator: 



(8) 



Here X[Si,Ti]it) is the characteristic function of the inter- 
val [Si,Ti], determining the duration of the pulse. The 
complex number C(j, depends on the state of the radiation 
as follows: 

cc, =2ttJ dkg*{k)fi{-k)S{uj{k) - oji) (9) 

and determines the Rabi frequencies of the field. 

In fact, the coherent vector is an eigenvector of 
the time rescaled annihilation operator aA,a)i 

(t) (defined 

by ^) with eigenvalue 



J2 Jdr^J dk5*(k)/,(k)e^(- 



t)((^(k)-(^,)/^ 



xe 



From that, using the limit (see |g for details) 

lim J_e»(^-*)('^(k)-'^i)/A'giT(t^,-..„)/A= 

A^O A2 



one obtains the limiting eigenvalue 

n ^™ 

y dT6{t-T) J dk5*(k)/„(k)J^ 



X[S,.T,](t)Cc. 



We will consider the case Si — S , Ti ~ T for all I, which 
means that all lasers are applied during the same time 
interval [S,T]. 

Let X be an observable of the atom and denote by 

{X)t = 

the average, over the limiting state of the radiation, of its 
time evolution. Taking the average of both sides of the 
quantum Langevin equation Q and using property Q 
one gets the following equation for the averaged observ- 
able: 



d{X)t 
dt 



mx))t 



+iXlS,T]{t) Y,{c^L^{X) - clL+{X))t. 



(10) 



This master equation for the case with presence of exter- 
nal field is one of the main results of the present paper. 
In the next section we will rewrite it in the equivalent 
form as an equation for the reduced density matrix. 



IV. THE QUANTUM MASTER EQUATION 

The dynamics of the atom can be described either in 
terms of observables, or equivalently, in terms of density 
matrices. The equation for the reduced density matrix 
can be obtained using the relation 

Tr(/5(t)X) = Tr(p(0)(X},) 

and equation HlOfl . It is a quantum master equation of 
the form: 



dm 

dt 



(11) 



The time dependent generator of this equation is given 
by the sum of its dissipative and Hamiltonian parts: 



£t(p) = -Cdiss(/5) - i[Hcs{t),p]. 



(12) 



As already stated in the Introduction, the dissipative part 
is the same as in the case the reservoir in the vacuum 
state: 



£diss(p) 



Re 7^ 2D^pD^ - pD^D^ - D^D^p 



(13) 

The Hamiltonian part is determined by the effective time 
dependent Hamiltonian 



H,s{t) = ^ Im7^i?+i?^ + X[s,T]{t){c*^D^ + c^D+) 



(14) 
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The first term in the brackets is the standard term which 
appears in the stocliastic hmit wlien tlie reservoir is in 
the vacuum state. It commutes with the free atomic 
Hamiltonian Hp^. The other terms are new and their 
presence leads to the important consequence that the ef- 
fective Hamiltonian at time t G [S, T] does not commute 
with the free atom Hamiltonian 



(15) 



Therefore the states of the atom driven by a long time 
laser pulse, which are approximated by stationary states 
of the master equation obtained from 1)11(1 replacing in 
the effective Hamiltonian H14|l X[S,T]{t) by identity, will 
not be diagonal with respect to Hp^. Our goal is to inves- 
tigate the structure of these states. In generic situation 
the interaction with the laser field drives the atom to a 
(dynamical) equilibrium state at an exponential rate. In 
this sense the interaction prepares a state of the atom. 
Thus our problem can be reformulated as follows: What 
kind of states can we prepare applying a laser pulse? For 
a long laser pulse these states are described by the time 
independent generator £ obtained from UlUI) . by setting 
S — > — oo, T — > -fcx), i.e. the characteristic function 
X[s,T] in l(Tn|) replaced by the identity. 

In the limit S — > — oo, T — > -|-oo the effective Hamilto- 
nian becomes time independent 



In this case the stationary state pst is a solution of the 
following equation 

C{pst) = Cdiss[pst) ~ i[HcS, Pst] = 0. (16) 

In the following sections we solve explicitly this equation 
for two and three-level atoms. 



The first term in the interaction Hamiltonian describes 
the process of emission of a photon due to transition of 
the atom from the excited to the ground state. The sec- 
ond term describes the opposite process, that is the ab- 
sorbtion of a photon by the atom and transition to the 
excited state. 

The real part 7 = Re of the susceptibility (|3J| deter- 
mines the inverse lifetime of the atom without interaction 
with external field and the imaginary part Alo = Inyy^ 
determines energy shift (cf. 0, sect. 5.5). The dissipa- 
tive part of the generator and the effective Hamiltonian 
for this system have the forms 

£diss(/5) = 7(2o-_po-+ - p(T+cr_ - o-+(T_p), 
HcS = Acj|l)(l| + f7*cr_ + ncr+, 

where J7 = c^^ is the complex Rabi frequency of the laser 
field. Equation (|15|l for a stationary state pst of the two- 
level atom equivalent to the following equations on the 
matrix elements p„in ~ ('7i|/3st|n): 

7Pii = Im(ripoi) 
Ji^Pw = i^ipii - Poo)- 

The general properties trp = 1, /5+ = p of the density 
matrix lead to the relations poo + Pii = l, Poi — P*o^ 
Pii — Pii- Using these relations one finds that, with the 
notation a = —ifl/j^, the stationary state is 



Pst 



Pii Pio 



Poi Poo 



1 



l + 2|a| 



\a\ a 



a* l+a 



(17) 



This state is pure only if a = 0. For very intense illumi- 
nation (|il|/|7t^| ^ 1) the quantity \a\ becomes very large 
and the atom becomes saturated with equal probabilities 
in the upper and lower levels, so that poo = Pii — 1/2, 
Poi = 0- For very low illumination, the stationary state 
is the ground state. 



V. STATIONARY STATES FOR A TWO-LEVEL 
ATOM 

Let us consider a two-level atom under the rotating 
wave approximation. 

The free Hamiltonian of the atom is 

i/A = eo|0)(0|+£i|l)(l|, 

where Eq and £1 are the energies of the ground and the 
excited states, so that the Bohr frequency a; = ei — Eq > 
is positive. The ground state is denoted by |0) and 
the excited state by |1). The transitions between atomic 
levels are described by the operators 

D = \0)(l\ = a^, D+ = \1){0\=<J+. 

The interaction Hamiltonian has the form 

Hint = i{<7- a+{g) - (T+ (g) a{g)). 



VI. STATIONARY STATES FOR A 
THREE-LEVEL ATOM 

The free Hamiltonian of a three-level atom is 

i7A = eo|0)(0|+£i|l)(l|+e2|2)(2|, 

where |0), |1) and |2) denote the ground state, the first 
and second excited states. For a non degenerate atom 
the energies satisfy the inequalities £2 > ei > eo- There 
are three positive Bohr frequencies uji = ei — Eq, '^2 — 
£2 — £1, ^3 = £2 — £0, which correspond to three possible 
transitions. The transitions between the atomic levels 
due to interaction with radiation are described by the 
operator 

Z?-dt|0)(l|+d;|l)(2|+dS|0)(2|, 

where dj are some complex numbers. In particular, the 
contribution of the first term in the sum to the interaction 
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Hamiltonian corresponds to emission of a photon by the 
atom and transition from the first excited level |1) to the 
ground state |0). 

For the three-level atom the three possible transitions 
are described by the operators 

Di := D^, = dl|0)(l|, D2 := D^, = d;|l)(2|, 

dl\0){2\. (18) 



This system of equations, together with the conditions 
tr p = 1, p+ = p, can be easily solved. With the notations 



Hence, with the notations 

I3 lu^,\dj? = cij+ibj, 



(19) 



(ttj and bj are the real and imaginary parts of 7^) the 
dissipative part of the generator, in the generic case, can 
be written as 



£diss(p) = 2 (aipn + 03^22) |0)(0| + 2a2P22|l)(l| 
- ai(/5|l)(l| + |l)(l|p) - ia2 + a3)(/5|2)(2| + |2)(2|/5), 

where pmn — (™|/5|?^) are the matrix elements of the 
density matrix. Introducing the Rabi frequencies of the 
laser field 

^j-=c^,dj, (20) 

the effective Hamiltonian can be written as 

i?cff = foi|l)(l| + (fo2 + 63)|2)(2| 
+Oi|l)(0|+172|2)(l| + r!3|2)(0| 
+flt|0)(l|+f^^|l)(2| + f^;|0)(2|. 

The matrix elements of the stationary state for arbitrary 
Rabi frequencies ilj satisfy the following system of equa- 
tions: 

aipii -I- 03^22 = Im (f^iPoi + ^3P02) 

(02 + 03)^22 Im (i^2Pl2 + f^3P02) 

(72 + 73)^02 + i^3{P22 - Poo) - i^2P01 + i^tPi2 = 

7iPoi + i^lipii - Poo) - 'i^2PQ2 + 'i^*3P21 = 

(71 + 72 + 73)^12 + i^*2{p22 - Pii) + i^iPm = i^lPw- 

Let us consider the case when the radiation, describ- 
ing laser pulse, is in a coherent state which drives only 
transitions from the ground state to the highest excited 
state, i.e., a laser field with frequency near W3. In this 
case rJi = ri2 = 0, r23 = 7^ is the only nonzero Rabi 
frequency and the system of equations for the matrix el- 
ements of the stationary state becomes simpler: 

aipii -I- a3P22 = Im {^Pm) 

(02 + 03)^22 = Im(rjpo2) 

(72 + 73)P02 + i^*{p22 - Poo) = 

7iPoi + iVL*p2i = 

(71 + 72 + 73)Pi2 - i^^*Pio = 0. 



-iVL 



02 



72 + 73 Ol 

the solution, which is the stationary state, has the form 

/ P22 P21 P20 \ 

Pst = P12 Pii Pio 

V P02 Poi Poo / 
/|ap a \ 

r|a|2 . (21) 

\ a* l-f |a|2 / 



1 



l + (2 + r)|a|2 



The obtained density matrix pst is the stationary state of 
the three-level atom in the laser field driving only tran- 
sitions from the ground state to the highest energy state 
of the atom. 

If the ratio r is very large, i.e. 02 3> ai, then the sta- 
tionary state describes the population inversion between 
the first two excited atomic levels. In this case the popu- 
lation of the first excited state of the atom is greater than 
that of the ground state. To find the physical conditions 
for this effect let us note that the quantity 



MiP / dk|5(k)p<S(c^(k)-c^i) 



determines the inverse lifetime of the atom at the level 
|1) in the free space. That is t ~ 1/ai. Hence, if the 
lifetime r of the atom in the first excited state is long 
enough, then we obtain the population inversion between 
the ground and the first excited level. 

As an example, consider a three-level atom in A config- 
uration, so that the transitions between the ground state 
|0) and the first excited state |1) are forbidden. This 
means that di = 0. In this case ai — 0, the lifetime of 
the first excited level is infinite, and the steady state is 
|1)(1|, i.e., in this case all the population is concentrated 
on the first excited state. 



VII. THREE-LEVEL LAMBDA-ATOM 

Let us now consider a three-level atom in A configu- 
ration in the general coherent field, that is both Rabi 
frequencies VI2 and Vl^ (|20|l are different from zero. 

The system of dynamical equations for matrix elements 
of the density matrix is (matrix elements are time depen- 
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dent, pnrn = Pnm{t), the dot denotes time derivative) 

Poo = 2a3/922 - 2Ini(rj3po2) 

/O22 -2(a2 + a3)p22 + 2Im(f22Pi2 + ^13^02) 

P02 = -(72 + 73)P02 - i^l{p22 - Poo) + i^2P01 
POl = i^2P02 - i^3P21 

P12 = -{I2 + 73)P12 - «f^2(P22 - Pll) + i^lPlO- 

The stationary solution of this system of equations is 

/ P22 P21 P20 \ 

Pst = P12 Pll Pio 

\ P02 Poi Poo / 
/O \ 

|r22p -niQ.^ 
Vo -n^ni |f^3l' J 

An important property is that it is a pure state: pst 
where 



VW2 



r!*|l)-O3|0) 



(22) 



In particular, if the ground state is doubly degenerate 
and if the coupling does not distinguish between the two 
ground states, then ^2 = and the unique stationary 
state is the dark state (cf. 1] for a discussion). 



VIII. CONCLUSIONS 

In the present work we derive, starting from exact 
microscopic dynamics and using the stochastic limit 
method, the quantum master equation 1)11(1 for an ar- 
bitrary atom driven by a general laser field. As an ap- 
plication of this equation, we study the cases of two and 
three-level atoms. The explicit form of the stationary 
states of such atoms are found. It is shown that for a 
three-level lambda-atom in two laser fields, which drive 
both allowed transitions, the stationary state of the atom 
is a pure state given by a superposition of the two lowest 
energy levels. Varying the intensities of the lasers, any 
such a superposition can be obtained. Hence the two 
lower levels of the atom behave like a single qubit and 
any state of this qubit can be obtained in a stable way 
and in exponentially small time. 
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